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Abstract 

Mobile communication channels are often modeled as linear time-varying filters or, equivalently, as time-frequency integral 
' operators with finite support in time and frequency. Such a characterization inherently assumes the signals are narrowband and 

I may not be appropriate for wideband signals. In this paper time-scale characterizations are examined that are useful in wideband 

, time-varying channels, for which a time-scale integral operator is physically justifiable. A review of these time-frequency and 

■ time-scale characterizations is presented. Both the time-frequency and time-scale integral operators have a two-dimensional discrete 
^ ^ ' characterization which motivates the design of time-frequency or time-scale rake receivers. These receivers have taps for both time 
^ I and frequency (or time and scale) shifts of the transmitted signal. A general theory of these characterizations which generates, as 

, specific cases, the discrete time-frequency and time-scale models is presented here. The interpretation of these models, namely, that 

■ they can be seen to arise from processing assumptions on the transmit and receive waveforms is discussed. Out of this discussion 
' a third model arises: a frequency-scale continuous channel model with an associated discrete frequency-scale characterization. 

(N 



Index Terms 

Time-Frequency, Time-Scale, Frequency-Scale, Delay, Doppler, Dilation, Doubly Spread, Time- Varying, Canonical Channel 
Models, Rake Receiver, Wideband Regime 

^ : I. Introduction 

It is common to assume that a received communication signal is composed of superpositions of different versions of the 
transmitted signal. These different versions arise from reflections of the signal off of scatterers in the environment. In the 
time-scale channel model, each reflection is a delayed and time scaled copy of the transmitted signal. The delays arise from 
differing path lengths from transmitter to scatterer to receiver Relative motion of the transmitter, scatterers, or receiver causes 
■ time dilations/contractions of the transmitted waveform x{t). Thus, each reflection is of the form. 
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^ and the received signal is a summation of the reflections characterized by C{a,b), the wideband spreading function^ , 
^ ' y{t) ^ ff C{a,b)^=x i*—^] dadb. (2) 



, We call a time-scale channel a wideband channel when the wideband spreading function has finite support. Due to the 
$-H ' physical limitations of signal propagation, it is reasonable to expect that C{a,b) has finite support. The maximum possible rate 
of change in path length, which is constrained by the speeds of the objects in the environment, limits the support of C{a, b) 
to a narrow range around the a = 1 line. Causality and the propagation loss associated with increasing path length effectively 
limit the support of C{a, b) to a finite range in the b direction. The support in the a direction causes a spreading in scale of 
the transmitted signal, and the support in the b direction causes a spreading in time of the transmitted signal. Thus, channels 
described by (2) are often referred to as doubly spread channels. 

Many signals and signaling environments satisfy the narrowband condition, an assumption under which the time dilations 
or contractions are modeled as Doppler shifts. Under this assumption, each received reflection of the signal is assumed to be 
of the form, 

Xr.eit) ^ ^{t - r)e^^'''' (3) 

In the narrowband channel model, the received signal is a superposition of time delayed and frequency shifted copies of the 
input and the channel is characterized by the narrowband spreading function S{6,t), 

y{t) =11 S{9, T)x{t - T)e^^''''*dTde, (4) 
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where S{9, r) typically has finite support in 6 and t due to the physical hmitations of the channel. The span of this spreading 
in time and frequency has proven to be a crucial parameter in communication systems [1-3]. Regardless of whether the support 
constraint is satisfied or not, (4) is a time-frequency description of a general time-varying linear system, 

y{t) = J h{t,T)x{t~T)dT. (5) 

When S{6, r) has no support constraint, the transmitted waveform and environment need not satisfy the narrowband condition. 

Kailath's pioneering work in his 1959 Master's thesis [4] and the concomitant development of the rake channel model 
provided a mathematical framework for capturing the energy associated with multiple transmission paths between transmitter 
and receiver using a discretization of the channel model. This work was furthered in 1963 by Bello, who proposed a discrete 
time-frequency characterization of the time-varying channel [2]. In [5], Sayeed and Aazhang reinterpreted this characterization 
from a diversity viewpoint, and used this canonical time-frequency channel characterization which combines a discrete set of 
time delayed and frequency shifted versions of the input signal, 

n=Ok=-K ^ ^ 

where 

S{e, t):^ // S{e', r')sinc {{t-t')W) sine {{9 - 0') T) er^'''^''-'^'^'^ AO' At' (7) 



to define a delay-Doppler RAKE receiver, a two-dimensional extension of the classic rake receiver The delay-Doppler rake 
takes advantage of the inherent added channel diversity associated with time-varying narrowband channels [5]. While the 
narrowband assumption is satisfied in many wireless communication signal environments, many wireless systems are wideband 
due to the higher data rates and multiaccess techniques [6]. Thus we may expect, in light of differences in the narrowband 
and wideband models, some advantages to the development of a canonical time-scale channel characterization in wideband 
communication scenarios. Motivated by this, [7, 8] used the channel in (2) to derive a time-scale canonical channel model 

m,n ^ / 

where ao, 6o are related to channel and signal characteristics, and 

Cm n = [[ 'C(a, 6)sinc ( m — ^ ) sine ( ?i ^ ) dadb (9) 

J J V In ao / V / 

An identical formula has been derived independently in [9, 10]. There is a difference, however, in the physical meaning of the 

decomposition in (8) between [7, 8] and [9, 10]. We will discuss this difference in Section IV where we will also present our 

point of view on canonical channel models. For us, a canonical model will refer to a time-varying linear system applied to 

a particular class of transmit signals whose output is measured through a particular observation procedure. For example, the 

time-frequency canonical model derived in [5] is based on bandlimited transmit signals observed at the receiver over a finite 

observation horizon (i.e., a time-limited receiver). As we discuss below, the time-scale canonical model can be derived from 

bandlimited transmit signals being observed at a scale-limited receiver. Furthermore, the new third canonical frequency-scale 

channel model introduced in this paper can be derived from scale-limited transmit signals being received at a time-limited 

receiver. We elaborate on this point in Section IV. 

Based on the above interpretation, in Section IV-A we introduce a frequency-scale time-varying channel model of the form: 

y{t)= I [ p(tJ,a)e^27rc^*^a; I - I dadtj (10) 

which is equivalent to (5) for positive time supported input signals and positive time horizon receivers, as we show in Appendix 
C. The canonical channel model derived from (10) is 

Vit) = Vc™,„e^™/(^-^^)l[T„T.](0«o^'^(«?}i) (11) 



where 



and Cm,n are coefficients which depend on the span of the observation time horizon (T2 — Ti), the scale domain bandwidth, 
and frequency-scale spreading function (see Equation (81)). 

Each of the three doubly spread canonical channel models discussed above motivates the development of a different two- 
dimensional rake receiver. A delay-dilation rake receiver based on the canonical time-scale channel characterization [9-11] 
leverages the diversity in wideband signaling environments in the same way that the delay-Doppler rake leverages the diversity 
in narrowband signaling environments [5]. Such a channel model and receiver may be particularly useful for ultra- wideband 
signaling due to the extremely wide transmission signal bandwidth [12, 13]. 
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A. Outline of paper 

In Section II we review background material on continuous narrowband (time-frequency) and wideband (time-scale) channel 
characterizations and examine simple one-path delay-Doppler and one-path delay-dilation channels in the framework of these 
representations. We derive and discuss the mapping between time-frequency and time-scale kernel operators and note that 
there exist time-frequency channels with no corresponding time-scale channel. In Section III we develop a general technique 
for the generation of canonical channel models and demonstrate the application of the technique to time-frequency and time- 
scale kernel operators. In Section IV we discuss the interpretation and derivation of these canonical models from reasonable 
processing assumptions on the transmit and receive waveforms. In Section IV-A we propose a frequency-scale canonical channel 
characterization based on the translation operators in frequency and scale. We conclude and propose future work in Section V. 



II. Continuous Narrowband and Wideband Channel Characterizations 

In this section we review and discuss the time-frequency and time-scale channel models and examine some simple channels 
to gain some intuition concerning the characterizations. The time-frequency description is a general time-varying linear system 
characterization. However, in a slight abuse of nomenclature, we wiU refer to all channel characterizations which can be 
related to the channel described by S{6,t) via Fourier transforms and phase factors as narrowband channels. Specifically, 
in this section, we discuss twelve such equivalent characterizations which were first explored by Kailath [4], Zadeh [14], 
and Bello [2]. We call these "narrowband" characterizations because when S{d,T) has finite support, the characterization is 
typically used only in narrowband systems and is not appropriate for wideband signals. We will only discuss the support 
condition constraint on S{9, r) for the narrowband characterizations when relevant, and consider the more general case where 
there is no such constraint on the support of S{d, t). Similarly, we will refer to channel characterizations based on the time-scale 
kernel C{a,b) as wideband characterizations because they are typically used in a wideband setting [15]. 



A. Narrowband Characterizations 

In this section, we develop a general technique for the generation of canonical channel models and demonstrate the application 
of the technique to time-frequency and time-scale kernel operators. 

The linear time-varying channel is characterized by the time-varying impulse response h{t, r) which denotes the response 
of the channel at time t to an impulse at time t — t. The channel input-output relationship is thus. 



y(t) - J hit,T)x{t-T)dT 

Such notation is used in, for example, [5, 16-19]. 

Another possible notation for the time-varying impulse response is 



yit) 



fco(i, T)x{T)dT. 



(13) 



(14) 



with the interpretation that ko{t,T) is the response of the channel at time t to an impulse at time t. This is the formulation 
used in, for example, [20-22]. Bello [2] calls fco(t, r) a kernel system function and notes the obvious correspondence between 
the two representations, h{t,T) = fco(t,t — r). Bello [2] defines four equivalent representations of the time-varying channel 
represented by ko{t, r) that map the time or frequency representations of the input into the time or frequency representations 
of the output. We define these four kernel functions. 



(15) 



The kernel system functions can be transformed into one another using the Fourier transform. For example, the kernel function 
that maps the input time domain to the output time domain (fco(i, r)) and the kernel function that maps the input time domain 
to the output frequency domain (ki{d,T)) are Fourier transforms of one another with respect to the first argument. We can 
summarize the relationships among the kernel system functions as follows. 



y{t) = J ko{t,T)xiT)dT Y{6) = J fci(0,T)x(T)dr 
/ k2it,iy)Xiiy)diy ¥(9)= I k:i{e,u)X{u)dv 



ko{t,T) 



kii9,T) 



V — >r 



(16) 
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h{t,T) 


input delay spread function 


y[t) = 1 h{t,T)x{t - r)dT 


s(e,T) 


delay-Doppler spreading function 


5(6»,t) = 1 h{t, r)e-J'2'^**dt 


T{t,u) 


time-varying transfer function 


T{t,v)= 1 r)e-J'2'^^''dr 


H{e,u) 


output Doppler spread function 


H(e,u)= II fc(t,T)e--'2'^(**+^'')didT 
,/,/ 


G(e,v) 


input Doppler spread function 


Y{9)= / G{e,v)X(e - iy)du 


V(t,v) 


Doppler-delay spreading function 




M{e,T) 


frequency dependent modulation function 


A/(6»,r)= 1 G{9,u)e^'^'"'^diy 


a{t,r) 


output delay spread function 


g{t,T) = II G(6»,iy)eJ'2'^('^*+''^'dediy 



TABLE I 

Eight system functions characterizing the linear time-varying channel, their function names from Bello [2], and their 

associated input-output relationship or definition. 



That is, 



(17) 

k2{t,v) = / /S3(6i,i/)e^27retj^ k3{0,iy) = j ki{0,T)e^^''^''dT 



The direction of the Fourier transform between fco and k2 (and also between fci and k^) is opposite to convention; We take 
the Fourier transform with respect to a "frequency" variable {i') and replace it with a "time" variable (r). This is necessary to 
be consistent with the kernel functions as defined in (15). 

Bello [2] provides the following useful interpretation of the kernel system functions, 

• The response to input 6{t — to) is time function fco(t,to) with spectrum ki{6,to). 



The response to input e-' " is time function k2{t,Oo) with spectrum k^{9,9[ 



0) 



and also notes, by simple inspection of (15), that fco and fca are time-frequency duals of one another, as are fci and fc2- 

Despite the simple input-output interpretations, the kernel system functions often lack intuitive physical interpretations [4]. 
For this reason, Bello [2] and Kailath [23] examined eight other system function characterizing the linear time-varying channel. 
These eight system functions are (13); its time-frequency dual. 



Y{e) = J G{e,v)x{e~v)di 



(18) 



the three functions obtained by taking the Fourier transform of h{t, r) with respect to t, r, and both t and r; and the three 
functions obtained by taking the inverse Fourier transform of G{9, v) with respect to 9, v, and both 9 and v. These eight 
functions are listed in Table 1. In the current literature, t) is usually referred to as the time-varying impulse response, 
(e.g., [5, 16-19]) and the delay-Doppler spreading function, S{9,t), is known simply as the spreading function (e.g., [5, 16- 
19, 21]). Unfortunately, fco(t, r) is also commonly referred to as the time-varying impulse response (e.g., [20, 21]). We will 
refer to ko{t,T) as the time-varying impulse response kernel to avoid confusion. 

The relationships among the eight functions via duality and the Fourier transform are summarized in the following diagram. 
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Duality is represented by a dotted line. 



h{t, t) 



T 

J T — >V 



M{0,t) 



T{t,u) 



(19) 



We can derive the following input-output relationships, 

y{t)^ J h{t,T)x{t^T)dT y{t)= jj S{e,T)ei'^'''>'x{t-T)dedT 

y{t)^ j T{t,u)e^'^'"'^X{v)Av Y{e)^ j H{e ~ v,v)X{u)<lv 

and 

Y{e)^ I G{e,v)x{e ~u)du Y{e)= u v[t,v)e-'^'^''^'^x{e-v)dtdv 

Y{e)^ J M{9,T)e-^^'''''^x{T)dT y{t) = J g{t - T,T)x{T)dT 

We can relate the eight system functions to the four kernel system functions as follows, 

ko{t,T) = h{t,t-T) ^g{t-T,T) 



(20) 



V{t, 6l)e^'2'^(*-^)(''+'')dTd6' 



(21) 



(22) 
(23) 

(24) 
(25) 



S{9,t) and ¥{1,1/) are distinctive in that their input-output characterizations and relations to the kernel system functions 
involve double integrals. In fact, it is the double integral formulation involving S{6,t) in (20) with the interpretation that 
the output is a superposition of time-delayed and Doppler-shifted copies of the input that makes S{d, r) an extremely useful 
characterization. For completeness, we note the inverse relations, 

S{9,t)= JJki{iy,t)e^^''^'+^^'^''-^Uiydt (26) 

Vit,iy)= JJk2{T,9)e^^''^*~^^'^'^+''UTd9 (27) 
and note the following relationship between the dual characterizations h{t, t) and G{9, v), 

h{t,T) = jj G{9,v)e^^''^'e-^^'''^'-^^'^^-''\\9dv (28) 

G{e,v) = jj /i(t,T)e-^^'^*''e^27r(e-,.)(t-r)j^j^_ (29) 

Although less commonly used in the literature, fc3((?, v) plays a pivotal role in understanding the narrowband and wideband 
characterizations [7]. We note the mapping between k^ and 5, 



k:i{e,v) = j S{9 ~v,T)e-^'^''^''dT 
S{9,t) = I fc3(6' + zy,i/)e^2^^''dz^. 



(30) 
(31) 
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Fig. 1. k(i{t,T) (left) and hit.r) (right) for the delay-by-ro channel. 



Fig. 2. S{9,t) for one-path channel with delay to and Doppler shift Bq (left); C{a, b) for one-path channel with delay feg ™d time dilation ao (right). 



which can be derived directly from the input-output channel characterizations. In the kernel system formulation of the channel, 
the outputs could be simply expressed in term of the kernel functions for inputs that were impulses in time and frequency. For 
the above characterizations, these relations are: 

• The response to S{t — to) is h{t, t — to) with spectrum M{d, to )e 

. The response to e^^'^'^o* is T{t, 9o)e^^'"*'^° with spectrum H{9 - Oo, 

For clarity, we display just the front face of the cube in (19), which details the Fourier transform relationships among the 
four most commonly used system functions. 

h{t,T) S{9,t) 



%). ' 



T{t,v) 



(32) 



H(6,v) 



In order to get some intuition concerning the channel characterization functions (both the four kernel functions (15) and 
the additional eight characterizations Usted in Table I, we examine simple channels. Consider the time-invariant channel that 
consists of a pure delay. 

X{t - To) 



x{t) 



channel 



In the case of the time-varying impulse response kernel, this channel is represented by ko{t,T) = (5(t — r — tq). In the case of 
the time-varying impulse response, this channel is represented by h{t, t) = (5(t — tq). Plots of these two functions are displayed 
in Figure 1 . One useful attribute of a system function is for a visual inspection of the function to readily reveal some physical 
properties of the channel. In the case of Figure 1, we see that, for ko{t,T), a diagonal delta function line crossing through 
(0, — To) and (to,0) arises from a delay of tq. For h{t,T), a delay of tq corresponds to a horizontal delta function line tq 
from the origin. A channel with several reflections (i.e., several different delays), would thus correspond to a system function 
with several parallel delta function lines. When the channel involves both a simple delay and a Doppler shift, the simple delta 
function lines for both ko{t,T) and h{t,T) are modulated by the Doppler shift. Table II displays the twelve system functions 
for the delay and delay-Doppler channels. The system function with the simplest form is S{9, r) which is the product of delta 
functions. From this, we interpret a region of localized energy in S{9, r) centered at {9o,to) as arising from an echo path with 
delay tq and Doppler shift 9o', see Figure 2. 
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y{t) = x{t - to) 


y{t) = x{t- To)e3'^^""^ 


koit,T) 


S(t - T - to) 


5{t-T- ro)eJ^'^'^o' 


ki{e,T) 


g-j27r(r + To)6l 


g-j2,r(T + To)(9-So) 


k2{t,y) 




gj2irt(i/ + flo)g-j2irToi' 






5{e~u~ 0o)e-J^'^^o'' 


h{t,T) 


S{t — To) 


5(r - ro)eJ^'"'''o 


s[e,T) 


5{t - To)S{e) 


6{t - Tn)5{e - Oo) 


T{t,u) 




gj2ir«otg-j27rToiy 


H{e,u) 


e-^'^^^o^sie) 


g-j2,rro^5(g _ Qjj) 


G{6,u) 




g-j2^To(H-tfn)5(,^_e),j) 


V{t,v) 


5{t - To)5(u) 


^j27VToOoS(t ~ To)5{u - 0o) 


AI{0, t) 


g-j27rTofl 


g-j27rro(e-So)gi27rr9o 


git,T) 


5{t - to) 


^f27VTaUaS(t - To)eJ^'"^''o 



TABLE II 

Time-frequency CHARACTERIZATION FUNCTIONS for the one-path delay and one-path delay-Doppler channels. S{9, t) has a very 

SIMPLE FORM FOR THE ONE-PATH DELAY-DOPPLER CHANNEL. 





y(t) = x{t - bo) 


y(t) = ^xC-^) 


C(a,b) 


5(a - l)5{b - bo) 


5{a - ao)S{b - bo) 


cc'\a,e) 


5{a- l)e-r^^i'<>» 


5(a-ao)e-J2^r6o« 



TABLE III 

Time-scale characterization functions for the one-path delay and one-path delay-dilation channels . 



B. Wideband Characterizations 

Starting from the wideband channel characterization, 



\a\ 

we derive the frequency domain to frequency domain mapping. 



y{t)= II C{a,b)^=x(- — - ] dadb. (33) 



Y{e) = III £{a,b)^=x(*-^^]e-'^^*'^dadbdt (34a) 



6) v/Hx (t') e-^'2'^(''*'+'')^dadWt' (34b) 
JJ £{a, b) a/HX {aO) e-^^'Tbedadfo (34c) 

£(2)(a,0)= f C{a,b)e-^^^''''db, (35) 



Y{e)^ d'^'>{a,e)y^\X{ae)da (36) 



and defining. 



we obtain 



Table III displays the wideband characterization functions for the one-path delay and one-path delay-dilation channels. In 
the narrowband case, S{9, r) is the product of delta functions for the one-path delay-Doppler channel; In the wideband case, 
the one-path delay-dilation channel is the product of delta functions. We interpret a region of concentrated energy in C{a, b) 
centered at {qq, bo) as arising from an echo path with delay bo and dilation parameter gq. 

C. Narrowband and Wideband Correspondence 

In this section we briefly examine the correspondence between the narrowband and wideband channel models. More 
specifically, we wish to link the narrowband channel model characterized by the dozen system functions discussed above, 
one of which was described by the time-frequency integral operator, 

Afsx{t):^ (( S{e,T)x{t~T)e^^''^^dTde (37) 
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to the wideband channel description embodied in the time-scale integral operator, 

6' 



C{a, b) — 7=x 



dadb. 



(38) 



We are interested in the mapping between S and C for Ms = Wc- The approach taken here differs from the traditional 
interpretation of the narrowband characterization as an approximation of the wideband characterization when applied to 
narrowband signals. This approximation is discussed in detail in, for example, [15, 24-29]. We do not consider the narrowband 
description of the channel as an approximation of the wideband channel, but rather look at the two descriptions without 
constraining the properties of the input signal. 

We first establish the relation from wideband to narrowband, showing that for every time-scale kernel, there exists a 
corresponding time-frequency kernel. Starting from (38), we have 



yit) = 



C{a,b)- 



1 



dadb 



and therefore. 



hit,' 



t - b 



\a\C{a,t — aT)da I x{T)dT, 



a\C{a, t — aT)da. 



(39a) 
(39b) 

(40) 



Returning to the mapping from C{a, b) to the narrowband characterizations, starting from (40), the remaining system functions 
can be related to £(a, 6) as follows, 

h{t,T)= / y^\C{a, {1 - a)t + aT)da (41) 



(42) 



and, taking the Fourier transform of (41) with respect to t, we obtain. 



S{9,t) = 



\a\C{a, (1 - a)t + aT)e-^^'''"dadt 



Using (40), it is possible to relate £ to all twelve narrowband representations [7]. 
It is also possible to express £{a,b) in terms of S{6,t), 



C{a,h) 



\0\ 



S{{1 - a)e,T)e^'^'''>^^-''^'^dedT 



(43) 



although the mapping relies on the assumption that the input signal has no DC component; see [7] for a discussion of this 
mapping. We can observe from (36) that, in the wideband model, the DC input component can only affect the DC output 
component. Intuitively, it is clear that rescaling the time axis and shifting in time a DC signal does not have any effect, and 
all the time-scale channel can do is amplify or attenuate the DC component of a signal. This is not the case in the narrowband 
model. For example, from k^{6,v) in (15) it is clear that the DC input signal component can affect any output frequency 
component. Therefore, there are time-frequency characterizations which have no corresponding time-scale representation. 

We look to some simple channel models and examine the mappings between C and S. We first consider the wideband 
(delay-dilation) single path channel, 

C{a,b)^5{a~ao)5{b^bn). (44) 



It follows from (42) that. 



- — -e 



|l-ao| 

6{e)5{T 



bo) 



ao = 1 



(45) 



and, substituting this into (4) we obtain y{t) — Xao^boit), as expected. 

We can derive the time-varying impulse response characterization h{t, t) for the wideband single path channel. 



et 



de 



|l-aol 

^\a^\ ^ fbo-aQT - {1- ao)t 



|1 - flol V 1 - ao 

\/\a^\5{ba - aor - (1 - ao)i) 



(46a) 

(46b) 

(46c) 
(46d) 
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one-path delay only 


one-path delay-Doppler, 7^ 


one-path delay-dilation, ao ^ 1 


S(e,r) 


5(e)5iT - to) 


s(e - eo)SiT - to) 


iY ' e i-^o 


C{a,b) 


5{a - l)S{b - to) 




<5(a — ao)(5(t> — 5o) 


h{t,T) 


<5(r-t,j) 


<5(r-ro)e-''2'^*^o 


^y\ao\5{{l - ao)t + aoT ~ bo) 




(5(6 - !/)e--'2irtoi. 


(5(e - J/ - 6»())e-327rTo!' 


V^|ao|e-J2'^«''o5(,/ - aod) 



TABLE IV 

Time-frequency and time-scale characterizations for the one-path delay-Doppler and one-path delay-dilation channels. 



which is also valid when ao = 1. We can compare this result to that of the single narrowband path (delay by rg, Doppler shift 
by 9q) channel, h{t, r) = 5{t — To)eJ^'"*^°. The wideband path gives rise to a delta function line with slope """^^ intersecting 
the r-axis at bo/ao; The narrowband path gives rise to a modulated delta function line parallel to the t-axis intersecting the 
T-axis at tq. 

We now turn to the expression of the narrowband (delay-Doppler) single path in the wideband model: 

S{0,t)=S{9-0o)S{t~to) (47) 
If we ignore the difficulties arising from the instabilities on the a = 1 line [7], it follows from (43) that 

-^la, ^) = — T=F-; ;^e^ (48) 

V|a|(l - ay 

and, substituting this into (2), we indeed obtain, y{t) = Xrafioit)- 

The various channel characterizations for the simple one-path models (including the time-invariant one-path model) are 
displayed in Table IV. We note that the one-path delay-dilation channel requires infinite support in time-frequency (45) 
whereas it requires only point support in time-scale (44). On the other hand, the one-path delay-Doppler channel requires 
infinite support in time-scale (48) whereas it requires only point support in time-frequency. Thus, since we are interested in 
channels which have finite support in time-frequency or time-scale (as we will see in the next sections), the choice of channel 
model is crucial and must be appropriate to the signaling environment (i.e., narrowband or wideband). Examination of ^3(6', ly) 
for the one-path channels reveals that it is possible (up to a scaling constant) for the one-path delay-Doppler and the one-path 
delay-dilation channels to have the same effect on a narrowband signal (eg, Xiiy) — S{i' — lyo)) by setting tq = bo/ao and 



ag 



III. Discrete Canonical Channel Models 



In this section we develop a general technique for the generation of canonical channel models and demonstrate the application 
of the technique to time-frequency and time-scale kernel operators. 

A. The canonical rake receiver model 

We begin with the derivation of the canonical model associated with the standard rake receiver. The classic expression of 
the sampling theorem for a signal X{v) with support {—W/2,W/2) is 

n— — 00 V W / 

An alternative formulation of the sampling theorem [19] is obtained by defining g{t) ~ x{a — t), 
and thus, 

nx sin(7rl^(t-f)) 



y ,f^)!^^H!L!!_AizMZ (50) 



n— — oo 



Mapping {a.,t) {t,T), we obtain. 



v-^ / n \ sm I 



sin {-kW (t-^)) 



(52) 
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Following [19], substituting (52) into the time-varying impulse response channel characterization (13), we obtain 



yit) 



h{t, T)x{t — r)dT 



sin(^l^(r-^^^ 



ttWIt 



W I 



(53a) 



(53b) 



(53c) 



n=0 



where the approximation is made based on the assumption that the channel is causal and has finite multipath spread, T^. That 
is, h{t,T) = 0,Vr < 0,T > Tm- Under this assumption, the approximation (53c) corresponds to hn{t) for which the mainlobe 
of the sine function overlaps with the support of the time-varying impulse response. The tapped-delay line in (53c) forms the 
basis for the classic rake receiver, where /i„(t)'s are usually assumed to be independent of each other 



B. The canonical time-frequency model 

We now proceed to examine the time-frequency canonical channel model which was originally derived in [5]. Alternative, 
but similar models are explored in [30-32]. The path we take in this derivation is essentially the same as that in [5]. We look 
at only the (0, T) portion of the received waveform, that is, j/(i)l(o,T)(i)- Starting from (53b), we impose the (0, T) restriction 
and obtain 



n 
W 



2/Wl(o,T)(t) = ^'(^ 

n— — oo 

Now we expand the r)l(Q (t) term as a Fourier series. 



j Mt,r)l(o,T)Wsinc (W (t 



n 
W 



dr 



/i(t,r)l(o,T)W 



5Z 2^ 

— — oo 
5Z T 



J27rkt/T 



fc=-c 



/i(t',T)e-^'2-'^*'/Tdt' 
/.(t',r)l(o,T)(i')e"^''^'*'/^dt' 



^327rkt/T 



(54) 

(55a) 
(55b) 



S(e,r)Tsinc((|,-0)T)e-3"('=-^<')de 



which is valid for t e (0, T). 

Substituting (55b) into (54) we obtain. 



n=—oo k=~cxi 



T' W 



where. 



S{e', r')sinc ((t - r') W) sine {{6 - 6') T) e 



(56) 



(57) 



(56) is valid for that part of any bandlimited signal received during (0, T). 

Under the path scatterer interpretation we assume that the channel introduces a maximum delay spread of Tm and maximum 
Doppler spread of Bd, that is, S{0,t) has support in {—Bd,Bd) x (0,Tm). In the smoothed version of S{9,t) in (57), if we 
consider only the terms in (56) where the main lobe of the smoothing kernel (which has size (— 1/T, l/r)-by-(— 1/W)) 
overlaps with the support of S{9, r), we need only sum over n ~ Q, . . . , N where N = [W^Tm] and k = —K, . . . ,K where 
K = iTBdl ■ We thus obtain the canonical representation of the time-frequency channel model, 

Y, E .(,-^).<-'/-5 14 (58, 

n=0 k=~\TBd'] ^ ^ 

C. Restatement 

The double sum time-frequency channel formulation (56) was obtained by assuming, 

• the input signal is bandpass with bandwidth W, and 

• the output signal is analyzed only for t € (0,T). 
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With these assumptions in mind, we define the following two projection operators, 

PTx{t) -.^ l[„^T]{tMt) (59) 

and, 

Qwx{t) ■.^T-'{l[_w/2,w/2]{uj)Hx{t)}{Lo)}, (60) 
and using the following two operators, the translation operator, 

TrX{t) -.^ X{t-T), (61) 

and the modulation operator, 

M^x{t) ■= x{t)e^^'"'\ (62) 

we can rewrite (56) as, 

PtUsQw = V c,„,„Pt A/rri Qw (63) 

where the c,„_„ ~ "^(t' w) -^s '■^^ narrowband channel operator defined in (37). Restating the channel operator 
in this setting, we can ask what general properties of the operators allow us to express the channel as a double summation 
of transformed input waveforms. In this section, we determine properties of the operators that are sufficient conditions for 
the existence of such an expansion. Our goal is to develop an analogous time-scale canonical channel model. That is, in 
Section III-F we propose projections P and Q such that, 

PyVLQ = ^rn^nPDZniQ (64) 

m,n 

for some choice of dilation and translation spacing parameters (uq and Bq), where the Cm,n depend on £, and D is the dilation 
operator, 

Dax{t) := -^x (-) , (65) 



for the wideband channel operator defined in (38). 

D. Generalization 

For the statement of the general theorem, we require the following definition. 

Definition 1 (paired-up operators): P and U are paired-up operators with generator bq iff, 

1) P is an orthogonal projection in L^(R) 

2) U is unitary in L^(M) 

3) PU = UP 

4) 3eo € RanP s.t.{t/"'eo : m G Z} is an orthonormal basis for RanP 

Using two different pairs of paired-up operators, the following theorem gives a sufficient condition for the channel expansion. 
Theorem 1: If {P,U) and {Q,V) are both paired-up operators with generator elements cq and /q respectively, iJ is a 
bounded operator, and 3c,„ „ such that 

J2 (y^^'fo, U'-"'eo) = (HV'fo, C/'eo) , Vfc, /, (66) 

then, 

PHQ = c^.nPU^'V'^Q (67) 

The proof of this theorem can be found in Appendix A and a method for calculating the coefficients „ can be found in 
Appendix B. 

E. Revisiting time-frequency 

The example we have seen so far of the application of this theorem corresponds to the situation 

. (P,C/,eo) = (PT,M^,^l[o,T](i)) 
• (Q, V, /o) = {Qw,Ti. , VWsinciWt)) 

w 

for the operator H = Ms of the form, 

Hx{t) =11 S{e, T)e^'2'^''*x(i - T)dedT. (68) 



Modulation and translation operators are a natural fit with our channel description, Afs, which describes the channel as a 
(continuous) summation of time and frequency shifts of the input signal. In Appendix B we demonstrate the coefficient 
calculation procedure for these specific operators. The procedure correctly derives the result Cm,n — S where S is 

defined in (57). 
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F. Time-scale canonical model 

We now develop the time-scale canonical characterization. For other possible extensions to time-scale, see the approach 
in [33-35] using wavelet packet modulation. 

The Mellin transform (also known as the scale transform) of a signal x E L^{0, oo) is defined by 

mx{uj):^ e-^'2"'"''^*x(i)-^ (69) 
Jo yt 

which represents the composition of two unitary transformations 

x{t) ^ e*/2a;(e*) mx{uj). 

For more information on the Mellin transform and its use in time-frequency analysis we refer the reader to [36]. For the 
time-scale canonical characterization, we will require the projection operator in the Mellin transform domain 

i?„ := 9K-H[-a/2,a/2]97l (70) 

which acts on a function x G L^{Q, oo) as follows 

x{t) ^ Tlxiiu) l[_„/2,„/2](c^)9Jtx(t^) R^x{t) 

where a > defines the cut-off MelUn "frequency". Explicitly, this means 

.CO ^ 

Rax{t) = / — =sinc[a(ln t - In T)]x(r)dT, t > 0. (71) 
Jo vtT 

Using the characteristic function in the Mellin transform domain, 

roH = lr 1 ^iH, (72) 

[ 2 In ao ' 2 In a„ J 

leads to the scale generator 

\ -ir Trsinc f : <> 

7o(i) = S vT^T^xA (73) 

[ : t<0 

For further details on the Mellin transform domain and its generators, consult [37]. 

It can be shown that (P, eo) = (R i , Daa , 7o(0) paired-up, and thus for the time-scale model, we use the following 

2 In ao 

paired-up operators. 



iQ,V,fo) = {Qi,n„,-^smcil)) 



to decompose the wideband channel corresponding to the operator H = Wc of the form. 



Hx{t) = // C{a,b)^=x (- — - ] dadb (74) 



into a discrete double summation. 



Py^LQ = J2 ^-n,nPDZT^,Q. (75) 



7n,n 



In Appendix B we calculate the coefficients in the time-scale case. 



/ [ C{a, b)smc ( m — ) sine ( n ^ ) dadb, (76) 

J J V Inao/ V «oo/ 



and the canonical time-scale model is then^ 

yii) = }:^-[--n^ ■ (77) 

m.n Uq \ V / 

-It is also possible to reverse the order of application of T and D in (75). In such case, the generated model in fact is more like that in (74) in that the 
time-shifts are not scaled. 
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IV. Physical Interpretation of Canonical Models 

In [10] the canonical model (77) is obtained by using two sampling results: the classical Shannon sampling formula for 
bandlimited functions, and a similar sampling result for functions that have finite support in the Mellin transform domain. 
With the help of these two formulas, [10] obtained a decomposition of the received signal into a series of such as (77) 
where parameters ao,bo are directly related to transmit signal bandwidth and transmit signal Mellin domain bandwidth. The 
trouble with such a model is that there are no signals that are perfectly (Fourier) frequency bandlimited and Mellin transform 
bandlimited (similar to the classical result that there are no time-frequency bandlimited signals except for the trivial zero 
signal). One can argue that the transmit signal is essentially frequency bandlimited, as well as, Mellin domain bandlimited, and 
thus a decomposition of type (77) should hold approximately. Furthermore, for a practical application, the infinite series (77) 
is truncated to a finite number of terms consistent with the finite size of the wideband spreading function C. Thus, another 
approximation is introduced, so overall one might expect that not much is lost by the initial assumption of joint Fourier 
frequency - Mellin domain band limitedness. 

In contrast, the approach we took in [8] does not suffer from the shortcomings outlined above. This different approach uses 
all the three players: the sender, the channel, and the receiver. The sender prepares the transmit signal by tailoring some of its 
properties. That is, the signal is embedded into the range of an orthogonal projection Q (e.g. Q can be an ideal lowpass filter); 
The channel acts via the operator H (5); and the receiver observes the channel output but in the observation process applies its 
own projection operator P through measurement, e.g. P is a time cut-off operator Thus, the entire transmitter-channel-receiver 
chain is modeled by a "sandwich" of operators PHQ where P and Q are under the user's control, and H is the channel 
operator. To simplify notation, the transmit signal is assumed to lie already in the range of Q, and thus Q often disappears 
from formulae. 

Thus, our basic model for transmitter-to-receiver communication system contains three blocks (see Figure 3): 

1) A transmit signal shaper, which is mathematically translated into a projection Q\ this can be thought of as the last stage 
of a modulator which, for narrowband communication channels, is either a bandpass filter around the carrier frequency, 
or a lowpass filter when analysis is done in the base band; 

2) A physical channel, mathematically modeled by a linear time-varying system H\ as such it can be written as in (5); 

3) A received signal observation, which again is translated into another projection P at the receiver; typically for memoryless 
source and channels, this is a time cut-off operator, due to real-time operation constraints. 

By changing the transmitter shaping and receiver observation projections, we obtain the different canonical representations. 
With this interpretation in mind, we can revisit previous models. 



X 



Q 



Transmitter 



H 

Physical Channel 



P 



Receiver 



Fig. 3. Our basic model for a communication channel. 



The standard rake receiver uses a channel model of type: 

y{t) ^Y.^n{t)x (i--^) 

n 

which is obtained for P equal to the identity operator (i.e. the entire channel output is available for processing) and for Q 
equal to the projection onto the space of frequency-bandlimited functions. As mentioned before, x is already assumed to be a 
frequency-bandlimited signal, thus x £ Ran Q. 

The time-frequency channel model of [5] uses the model (6) which is obtained when P is the time cut-off multiplication 
by l[o,T] ^nd Q is the ideal lowpass filter 

The time-scale channel model of [8] in (77) uses the ideal Mellin domain lowpass filter as P and the ideal lowpass filter as 
Q. In other words, the channel output is observed through a scale filter defined using the Mellin transform. We now present 
another canonical model in which the pair of projectors consists of the time cut-off 1[Ti,T2](^) for P and the ideal Mellin 
domain lowpass filter for Q. 

A. The Frequency-Scale Canonical Model 

We now consider a frequency-scale canonical channel characterization based on the translation operators in frequency and 
scale. In the frequency-scale model, we restrict ourselves to x{t) defined for t > and use for the transmitter projection Q the 
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TABLE V 

Summary of canonical models of the form PHQ = J]^ „ Cm,nPU"^V"Q. 



Mellin domain band limiter R i . Thus, the transmitter transmits scale-limited waveforms. We use for the receiver projection 

2 In aQ 

P simply a time cut-off 

P[Ti.T,]x[t) :=l[T„T2](t)a:(0 (78) 
where T2 > Ti > define the receiver observation time horizon. The overall chain of operators then decouples into the series 

PHQ -^Y.^"^^-PlT.-T.]Mr/^T,-Ti)D:^R^- (79) 

For this model the following theorem gives a decomposition into a series of dilated and frequency shifted versions of the input 
signal. 

Theorem 2 (The Canonical Frequency-Scale Channel Model): Assume a time-varying channel H defined by (5). Then for 
any signal x that is Mellin domain bandlimited to I—ttt — , — 1, i-C- x £ Ran Q, 

y{t) := Hx{t) = Y: c„,„e^^™"*4j. f±) (80) 

for all Ti <t< T2, where n ^ t^^, 

Cm,n = i^e-^-"(^^+^=) p ^j^ p{^, a)e^--(^^+^^)sinc - m) sine - da) dc, (81) 

and p is computed in turn from h{t,T) through (113). 

The convergence in (80) is in the sense. The proof of Theorem 2 is included in Appendix D. In (80) we see that if we 
receive scale-limited waveforms over a finite time window, that we can decompose the time varying channel into a discrete 
representation involving a countable sum of weighted scale frequency shifts of the transmitted waveform. 



V. Summary and Future Work 

Table V summarizes the projection and translation operators used to generate the three discrete canonical channel model 
discussed in this paper Each of the models can be thought of as sending a transmit waveform through a shaping transmission 
filter Q and then receiving the signal through a receiving filter P. The corresponding discrete channel models are presented 
in Table VI. We have presented here a general theory which generates these models based on assumptions on the transmitter 
and receiver characteristics: 

• The time-frequency model arises from: 

- frequency bandlimited transmit waveforms 

- put through a time-frequency (narrowband) channel 

- at a time-limited receiver 

• The time-scale model arises from: 

- frequency bandlimited transmit waveforms 

- put through a time-scale (wideband) channel 

- to a scale-limited receiver. 

• The frequency-scale model arises from: 

- scale-limited transmit waveforms 

- put through a frequency-scale channel 

- to a time-limited receiver. 

One of the many items for further study is the question of the physical interpretation of the frequency-scale model. In what 
settings can we envision a channel which imparts a limited range of frequency and scale shifts of an input signal? Perhaps a 
direct path only model of a wideband sonar signal reflecting off the undulating surface of the ocean with a moving transmitter 
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TABLE VI 
Summary of canonical models. 



would impart simultaneously a frequency shift (caused by the frequency of the ocean surface waves) and a scale shift (caused 
by the change in transmission path length during transmission). Indeed, one main topic of future research is to characterize 
the channel scenarios which lead to efficient representation in each of the three models. 

Further research topics include a full analysis of the two dimensional delay-dilation and Doppler-dilation rake receivers 
which arise from these canonical models, including an analysis as to which communication scenarios result in performance 
gains for the two-dimensional rake over conventional receivers. Also, we hope to generalize the information theoretic analysis 
to the delay-dilation and Doppler-dilation rake receivers similar to that which was done for delay-Doppler rake receiver 
in [5]. Similarly, it would be of interest to develop a canonical time-scale and frequency-scale multiantenna wideband channel 
model similar to that proposed in [38] for time-frequency channels. Also, [33-35] introduce wavelet-based channel models; A 
comparison of these models to the model derived in this work in Section III is a topic of future research. Finally, we ask, is 
there a corresponding underspread/overspread theory (see [22, 39]) for the time-scale and frequency-scale canonical models? 

Appendix 
A. Proof of main theorem 
Proof: First we expand PQ using the orthonormal basis and unitary properties of the paired-up operators, 

P-^(-,C/'"eo)C/'"eo (82) 



and 



we derive, 

PQx = "^{Qx,U"'eo)U"'eo (84a) 

m \ n I 

= 5](x,l/"7o)(^"/o,C/'"eo)C/™eo. (84c) 

rn.n 

We use this to determine, 

i'fjIcm.nC/"^"] gx - ^ c,„,„C/"PQy"x (85a) 

\m,n / m,n 

- Cra.nU^ j (^'^0' ^^^o) (^"^' ^'/o) ^'^0 ) (85b) 

m,n \ k,l I 

= c™,„("t^7o,C/'eo)(x,y-"l/'/o)[/"[/'^eo (85c) 

m.n.k.l 

= E (E (^"^"/o, C/-'-™eo) ] (x, y"/o) C/^eo (85d) 

where the commuting property of paired-up operators was used in (85a), (84c) was used in moving from (85a) to (85b), and 
the unitary property of V was used in moving from (85b) to (85c). Now, looking to the LHS of (67), we expand using the 



orthonormal basis and obtain, 

PHQx = {HQx, U'eo) U'eo 



Given H, we then compute, 
which we use to solve, 

for Cm,n- These Cm,n satisfy (67). 



= ^ {x, Wo) {HV^fo, U-'eo) U'eo 

s.u 

hu^, := {HV^fo,U'eo) 



B. Solving THE COEFFICIENT EQUATION 

We now discuss the form of the solution to (66). We define 



and define 



which allows us to express (66) as. 



*~"m^n • — '~"n, — m 

hu,s = 5Ic™."(^"+"/o,f/-'"'"eo) 
= ^(F"-"/o,C/^-"eo)£„,™ 



where 



k,l k,l 

Expressing h, a, and c in the Z-transform domain, 

^(zi,Z2) :-Em^i4«M -5]^f4(l^'/o,C/'eo) 



k.l 



we can write (91c) as, 
and solve for C 

In terms of Cm,n, this is. 



C{zi,Z2) :=Efc,;^i4cfcJ 

if = AC 



C(2l,Z2) 



A{zi, Z2) 
_^ ( H{zi, Z2) 



where 



Z 



(F(zi,Z2))„„ = / / e-^'2^^i™e-^'2^*^"F(e^'2^^Se^'2'^''=)d6li 
Ja Jq 
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and 



We can express (98) as a convolution of coefficients by defining 

^(gj2^ei J2^9.) .^ 1 (100) 

a„,„:^/ / e-J2''''i™e-J^''^^"i(eJ'2'^''Se^2^''^)d6iid6i2, (101) 
Jo Jo 

and we can obtain the „ using 

Coefficient calculation 

Thus, to calculate the coefficients Cm,n, 

1) calculate /ifej via (87), 

2) calculate a„i.„ via (89), 

3) use a,„,„ to obtain A(e^^''^S e-'^'^^^) via (93), 

4) use A(e-''2"''i,e^27re2) obtain a,„,„ vua (100) and (101), and 

5) use hk.i and „ to obtain Cm „ via (102). 

Here, we present the highlights of the coefficient calculation procedure for the time-frequency and time-scale canonical 
models. For more detailed steps, consult [7]. 



Example: time-frequency 



hkj = JJJ l[o,T](i)e^'"*''"*^sinc(Wi-fc-T^T)5(0,r)dMrdt (103) 
ara^n = ^J^ e-'^""^ smc{Wt - m)dt (104) 



For 9i,d2e [0,1], 



^.j2.e. „ f VWTe^^-™^ : 0^ G (O, i) 



dm.n = -^77^ I e-'^^'^'^mciWre^ + m) - J^^a_™,„d02 (106) 
\WT Jo yy 1 

c„ „ = jj S{9,T)e^''^'^''+"''>smc{Te + m)smc{n + WT)d9dT (107) 

which are precisely the coefficients in (58). 



Example: time-scale 



hu s = , // — i=C(a,b) ( f ^sinc | — z u ) sine ( s ) dt \ dadb (108) 



bolnao Jq y/t \bo J Vl^ao 
For 6*1, 6*2 G [^5' 5]' ill distributional sense. 



1 ft \ f\n\t\ 
sine m sine n ] dt (109) 



A{0,,92)^J^^bl^'''''^ f t-"^+i^^^oe^^^'^'dt (110) 





,'5 b ^ ^""0 p-j^TTdim -j2-!Te2n 

arn.n^V^ / 1 dOidO^ (111) 

£(a,5)sinc ( m — ] sine (n — ) dad5. (112) 



2 2 



Inao/ V ab, 



"0 
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C. The equivalence between (5) and (10) 

In this section we obtain the correspondence relations between the two forms (5) and (10) of a general time-varying linear 
system when input signal are supported on positive time domain, and the observation is restricted to a positive time horizon. 
Consider first the input-output relationship given by (5). For positive time supported input signals, the output is given by 



y{t) = / h{t,t-T)x{T)dT 

Jo 

We change the integration variable t — > ^, and since we have a positive time horizon, i.e. i > 0, we obtain: 

Now denote p(t, a) ~ ^-^/i [t, t — ^-j, and p{uj, a) its Fourier transform with respect to t. Then the inverse Fourier transform 
allows us to write 



yft) 



that is (10), where, explicitly. 



p{uj,a)^-^[ e-^^'"^Hh(t,t- -\ dt (113) 

a^/a Jq \ a J 

For the converse, assume the input-output relationship in given by (10). Then, performing the integration over uj first we 
obtain ^ 

y{t) = / p{t,a)—x ( - ) da. 
Jo va V«/ 

Next we need to change the integration variable a into t ~ t — - 



y{t) - / p[t, J ^—^x{t - r)dT 



which is exactly (5) with 
where lt>r(i) := l[r,oo)(0- 



h{t,T) = Uyrit)\l^ j e^^^'^'p[iO,j—] diO (114) 



D. Proof OF Theorem 2 

We follow the recipe proposed in Appendix B. The two sets of paired-up operators and generators are: (P = 1 [Ti .T2] j , cq {t) 
v^l[Ti,T2]W) and {Q ^ M-^l[__i_ _i_]m, DaoJoit) = thn;^sinc(j|^)lt>o(t)). First we need to compute /i^,; and 



flm.n- We have: 

foc ( poo I i*T2 



Next we compute 22) = E™ „ am,n2r^2 > at 21 = eP'''^^ Z2 = e^^^^^^ for 6*1 e [-\, i], 6I2 G [0, 1]. We obtain: 



V (fc'2 + m) Inao 



where no = no(6'2) is the only integer so that ^''+"" e [Ti, r2)- Then 



f2 

which has its Fourier expansion with coefficients a,„_„ given by 

rT2 



A(p32^ei „j2.e2) ^ ^ ^ / (g2+no)ln(^ j-2^e,^ in 



^/b^ / ' Vte^-''2"""*sinc (m+ ^\ dt 
Jti \ In ao / 



Then the coefficients c„i.„ that solve the equation a-kf ~h with = Cm,-n are given by 

c™,„ = [a * ;i)_„,„. = J_e-^"-^^(^^+^^) a)e^---(^^+^^)sinc - m) sine - da) d^ 

which is exactly (81). 
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